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We present preliminary results on extractions of the chiral LECs L\q and C%-/ and constraints on 
the excited pseudoscalar state 7r(1300) and 7r(1800) decay constants obtained from an analysis 
of lattice data for the flavor ud light quark V-A correlator. A comparison of the results for the 
correlator to the corresponding mildly-model-dependent continuum results (based primarily on 
experimental hadronic T decay data) is also given. 
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1. The V-A correlator 

We focus on the difference of flavor ud vector (V) and axial vector (A) current-current 2-point 
functions, n^J A , and their J = 0, 1 scalar components, ITyj A , defined in Minkowski space by 

n£/> 2 ) =ij d'xe^(0\T (4 /A (x)4* A (0)) |0> 

= (?V-*V v )n« (<2 2 ) + ?V<) A (<2 2 ) (l.i) 



Q 2 > 0, also determine the corresponding Euclidean 2-point functions 



V/A^ ) 1 1 V ±± V/A\ 



where, as usual, Q 2 = -q 2 . In what follows, we denote ATI^ = U ( y ] -U A J) . The U { y) A {Q 2 ), for 



u v/a(Q ) „..„, = (2 2 ^ V -G^G V ) n^ A (G 2 ) - G^G v n^ A (G 2 ), (1.2) 



Eucl 



making the Uyj A (Q 2 ) accessible from lattice simulations. Tl A ^ and U A ^ both have kinematic poles 
at Q 2 = while the / = + 1 sum does not. Since, beyond NLO in the chiral expansion, the 
pole residues involve at-present-unknown NNLO LECs, we focus on An(£ 2 ) = nj? +1) (G 2 ) - 
ni° +1 '(G 2 )> which satisfies an unsubtracted dispersion relation with only physical singularities. 
The corresponding spectral function, Ap(s), consists of a 5-f unction at s = m| and continuum be- 
ginning at 4m 2 . The "continuum part", AH(Q 2 ), of AYl(Q 2 ) results from removing the % pole: 
An(<2 2 ) = An(2 2 ) — 2f 2 1 (Q 2 + m 2 ). For s < m 2 , Ap (s) can be determined experimentally from 
hadronic % decay data [jl]]. Public versions are available from both ALEPH [Q] and OPAL 
with a yet-to-be-corrected problem affecting the covariance matrix of the former. Beyond s = m 2 , 
an alternate representation results from fitting a physically motivated model for duality violations 
(DVs) ^ to integrated versions of the data [Qj. The OPAL data and fitted DV model provide a dis- 
persive determination of An(<2 2 ) at spacelike Q 2 > 0, where it can also be measured on the lattice. 
The dispersive result is nominally quite precise, but has some (mild) model-dependence from the 
use of the DV model. It also involves X — > \%v z contributions in the V channel whose uncertainties 
may have been underestimated, given that the corresponding 4n branching fractions differ from 
expectations based on CVC and measured e + e~ — > A% cross-sections by much more than is typical 
for isospin-breaking corrections [^| 0|. The mildness of the model-dependence follows from the 
fact that, in the range Q 2 < (500 MeV) 2 expected to be of relevance to the determination of chiral 
LECs, DV contributions to Ap(s) account for a few to several % of AI1(<2 2 ), the precise values de- 
pending on the point in the spectrum at which one switches from data to the fitted DV model. The 
cross-check on the lattice results for All provided by the dispersive representation is also useful 
in light of the freedom to vary m u j iS on the lattice, which, in principle, provides access to NNLO 
chiral LECs currently unknown and/or difficult to extract reliably with continuum methods. 



2. Lattice data for ATl(Q 2 ) 

AI1(£2 2 ) has been determined for the fine 1 /a = 2.28 GeV, m n = 289, 345 and 394 MeV, and 
coarse 1 /a = 1.37 GeV, m % = 171 and 248 RBC/UKQCD DWF ensembles detailed in Refs. 
The latter provide an increased number of low-<2 2 points, improving the determination of the chiral 
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LECs. The values of f n and m n needed to convert An(<2 2 ) to AH(Q 2 ) are given in Refs. [[| |9]]. 
The simulation m s values are, in all cases, close, but not exactly equal to the physical value. As the 
SU (2) LECs correspond to physical m s , we analyze the data in the chiral SU (3) framework. The 
ensemble values required for this puipose are also given in Refs. The spectral function, 

Ap, of AIT, and hence also AIT itself, are 0(m9) in the chiral expansion. We thus expect the lattice 
AII(<2 2 ) to approach the physical results for sufficiently light ni(. We find that, within errors, for 
the low-<2 2 region of interest to us here, the lattice AYl(Q 2 ) agree well with one another, and with 
the continuum OPAL+DV model results, for all but the m n = 394 MeV case, as shown in Fig. [j] 
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Figure 1: Lattice and OPAL+DV model results for AIl(g 2 ) 



3. The chiral LECs 

At NLO in the chiral expansion, All is controlled by the single NLO LEC L r m (pL). Two 



previous lattice studies determined L r w {jx) by analyzing AIl^ 1 ' [10, 11] at NLO. With the lattice 
spacings available, the second-smallest non-zero Q 2 were found to be too large (~ (650 MeV) 2 and 
~ (460M<?V) 2 for Refs. [£(J and [|l|] respectively) to allow for a successful NLO fit. Final NLO 
analysis results were thus based on the single lowest Q 2 values, (320 MeV) 2 and ~ (230 MeV) 2 , re- 
spectively. The current analysis improves on the previous ones in a number of ways. First, the new 
coarser lattices allows access to an increased number of low-<2 2 points. Second, the statistics for 



the m n = 289 MeV ensemble considered previously [11] have now been doubled. Third, because 
the residue of the Q 2 = kinematic pole in An' 1 ' involves an unknown NNLO contribution, whose 
contribution, relative to that of the term involving the NLO constant U m , gets enhanced when one 
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goes to the low Q 2 desirable for extracting LECs, we switch to analyzing All rather than All' 1 ', the 
71 pole contribution which must be subtracted to obtain the former having both an exactly known 
residue and being farther from the region of the lattice data than is the kinematic pole in AIT 1 ). 
The NLO results for L[ (n), for pL = Ho = 0.77 GeV, obtained for each Q 2 < 0.25 GeV 2 , and all 
but the heaviest m % = 394 MeV ensemble, are shown in Fig. ^[ Also shown for comparison are the 
results of an NLO analysis of the OPAL+DV model over a similar range of Q 2 . 

Since the <2 2 -dependence of Hy(Q 2 ) is known to be poorly reproduced by the NLO represen- 
tation [12], one might be surprised by the relative stability of the results for different Q 2 . However, 
the missing intermediate p contribution believed responsible for the NLO Tly slope problem [|l2|] is 
encoded in the NNLO LEC C' 93 [|T|], and the contribution to Hy(Q 2 ) proportional to C' 93 is exactly 
cancelled by that proportional to C r 93 in the NNLO expression for II^ +1 \q 2 ). Nonetheless, the cen- 
tral value for the average slope of AII(<2 2 ) with respect to Q 2 is significantly larger than expected 
from the NLO expression, albeit at only the ~ 2a level. The structure of the full NNLO result, 
known from the results of Ref. fll3[], is very linear in Q 2 for the Q 2 considered here, so a significant 



portion of the NNLO contribution is easily removed by fitting the results of Fig. g to a linear form 
and using this to extrapolate to Q 2 = 0. The only NNLO contributions remaining to be removed 
are then those entering AIT(O). These involve two NNLO LEC combinations, Cgj — C[ 2 — C%o> 



which is not large-A^ suppressed, and C 62 —C[ 3 —C' 8l , which is [ |13| , |14| ]. The first combination has 
been estimated in Ref. [|l4|] using the results of previous continuum works. The second is currently 
unknown, and has only been loosely bounded, using rough large-A^-suppression arguments [pp. 



The resulting NNLO LEC combination assessments were used in Ref. [ 14 ] to obtain a continuum 
extraction of L' 10 and C' gl (the NNLO LEC expected to dominate the slope of All). The analysis 
employed the ALEPH data and was based on two additional, not explicitly tested, assumptions, 
namely (i) that the NNLO form will successfully represent AYl(Q 2 ) and (ii) that the V and A chan- 
nel DV spectral contributions which, being governed by the resonance structure in the channel in 
question, are expected to be channel-dependent, can be assumed to be approximately the same in 
form and hence combined into a single DV ansatz for the V-A difference. While the latter assump- 
tion is not borne out by the combined V and A channel fits of Refs. [Q], the contribution to AIT of 
the DV part of Ap is small in the low-<2 2 region. It is thus of interest to compare the results of our 



fit to those of this mildly model-dependent continuum analysis, which are [14] 



L\M = -0.0041 (4) NNL0 ; C^o) = 0.0049(2)^ LO GeV~ , (3.1) 



with the error dominated entirely by the uncertainty in the estimate for unknown large-A^-suppressed 
NNLO LEC combination. Results obtained using instead the OPAL data, and incorporating the re- 
sults of the DV model fits of Refs. [Q] in order to remove the second of the two assumptions noted 
above, are in extremely close agreement. An additional unknown systematic error, associated with 
the two additional assumptions noted above, of course also exists for the results of Eq. (3T). Our 
final goal is to perform a NNLO analysis of the lattice data, including a range of m q sufficient to put 
constraints on the currently unknown NNLO LECs (something likely to be feasible given the differ- 
ence of the results for AIT for m % = 394 MeV from those for the smaller m % seen in Fig. [I]), but this 
analysis has not yet been completed. The results of following the continuum estimates/arguments 
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of Ref. []T4|] for the unknown NNLO LEC combinations entering AI1(0) are 

L\M = -0.0038(4) /fl „ (4) NNL0 ; Cg 7 (jUo) = 0.0040(21) /fl „ (2) NNL0 GeV- 2 , (3.2) 

in good agreement with, though less precise than the NNLO-LEC-induced part of the error obtained 
explicitly in the continuum analysis of Ref. [14]. 
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Figure 2: NLO results for L' 10 (/Xo) from analyses of lattice and OPAL+DV model versions of ALT 



4. Constraints on the %' and n" Decay Constants 



Excited 1=1 pseudoscalar mesons, P, couple with strengths 2fpmp to the divergence of the 
flavor ud axial current. Their decay constants, fp, enter the conventional determination of m u + nid 



employing sum rules for the two-point function of this divergence [16] and are currently determined 
as part of the analysis. A similar internal determination is required in the extraction of m u + m s from 
sum rules for the two-point function of the divergence of the flavor us axial current. Finally, flavor- 
breaking sum rules used to determine | V„ v | from hadronic X decay [[l7|, or a combination of 
hadronic X decay and electroproduction cross-section data [|9|], encounter a problem with the very 
bad convergence of the OPE representation of J = contributions, necessitating the subtraction 
of the chirally suppressed, but not totally negligible, strange excited state scalar and pseudoscalar 
contributions to the differential X decay spectrum. The strange pseudoscalar subtraction relies on 
the excited K decay constants obtained in the sum rule analysis. 

The lattice data allows us to test the reliability of the sum rule determination of such de- 
cay constants, as follows. The quantity P(Q 2 ) = Q 2 AYl^°\Q 2 ) , which, for m u = is equal to 
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— <2 2 n^°'(<2 2 ), is free of kinematic singularities and satisfies a once-subtracted dispersion relation. 
Since H^\q 2 ) and the quantities m K , f n which determine the pion pole contribution to the dis- 
persive representation are all measurable on the lattice, the following rearranged version of this 
relation provides constraints on the continuum contribution, and hence on the excited state decay 
constants, f n i and f K », for each pair of Q 2 and subtraction point 

Spectral positivity ensures that the LHS provides an upper bound on the contributions from any sub- 
set of the full set of excited pseudoscalar states. In the narrow width approximation, this constraint 
represents a straight line in the f%i-f%n plane for each pair (Q 2 , Qq). The fact that the excited state 
decay constants scale as m 2 can be used to scale each such bound from the masses used in the sim- 
ulation down to physical m n . It turns out that, at present, only the high-statistics, l/a = 2.28 GeV, 
m % = 289 MeV ensemble provides data sufficiently accurate for this purpose. The envelope of 
the resulting set of constraint lines, scaled down to physical m % , is shown in Fig. |3| Also shown, 
for comparison, are the results obtained/used in Refs. [16]. These are obviously in good agree- 
ment with the lattice constraints, leaving room for small additional contributions to the dispersive 
representation from yet higher excited pseudoscalar resonances. 




Figure 3: Lattice constraints on the %' and n" decay constants 
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